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ABSTRACT 



We present a measurement of the rate of type la supernovae (SNe la) from 
the first of three seasons of data from the SDSS-II Supernova Survey. For 
this measurement, we inchide 17 SNe la at redshift z < 0.12. Assuming a 
fiat cosmology with Qm = 0.3 = 1 — Qa, we find a volumetric SN la rate of 
[2.93lo;j4(systematic)l:o;7'i'(statistical)] x 10~^ SNe Mpc"^ /i^g year"\ at a volume- 
weighted mean redshift of 0.09. This result is consistent with previous measure- 
ments of the SN la rate in a similar redshift range. The systematic errors are 
well controlled, resulting in the most precise measurement of the SN la rate in 
this redshift range. We use a maximum likehhood method to fit SN rate models 
to the SDSS-II Supernova Survey data in combination with other rate measure- 
ments, thereby constraining models for the redshift-e volution of the SN la rate. 
Fitting the combined data to a simple power-law evolution of the volumetric SN 
la rate, ry oc (1 + z)'^, we obtain a value of (3 = 1.5 ± 0.6, i.e. the SN la rate is 
determined to be an increasing function of redshift at the ~ 2.5a level. Fitting 
the results to a model in which the volumetric SN rate, ry = Ap{t) + Bp{t), 
where p{t) is the stellar mass density and p{t) is the star formation rate, we find 
A = (2.8 ± 1.2) X 10-^^ SNe Mq^ year-\ B = (9.3l|f) x 10"^ SNe M^K 



Subject headings: supernovae: general 
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Introduction 



Type la supernovae (SNe la) have gained increasing attention from astronomers, primar- 
ily due to their remarkable utility as cosmological distance indicators. There is now broad 
consensus that a type la supernova is the thermonuclear explosion of a Carbon-Oxygen white 
dwarf star that accretes mass f rom a binary companion until it reaches the Chandrasekhar 
mass limit (e.g. [Branch et al.l (119951 )). However, much remains to be learned about the 



physics of SNe la, and there is active debate about both the nature of the progenitor sys- 
tems and the details of the explosion mechanism. For example, the binary companion may 
be a main-sequence star, a giant or sub-giant, or a second white dwarf. The type of the 
companion star determines in p art the predic ted time delay between the formation of the 
binary system and the SN event (lGreggioll2005l ). The time delay can be constrained observa- 
tionally by comparing the SN la rate as a function of redshift to the star formation history 
(SFH) Jstroker et allbooi ICappellaro et allboOTh . 



The insight into the nature of the progenitor systems that SN la rate measurements pro- 
vide can also potentially strengthen the utility of SNe la as cosmological distance indicators. 
Although the strong correlati on between SN la peak lumin osity and light curve decline rate 

9931) . the physics underlying t his rela- 



was found purely empirically ( Pskovskii 1977 : IPhillips 



tion has been extensively studied (IHoflich et al. 



1995 



19961 iKasen fc WooslevI 120071 ). There 



is hope that improved physical understanding and modeling of SN la explosions, coupled 
with larger high-quality observational data sets, will lead to improved distance estimates 
from SNe la. As part of this program, deeper understanding of the nature of the progenitor 
systems can help narrow the range of initial conditions that need to be explored in carrying 
out the costly simulations of SN la explosions that in principle predict their photometric 
and spectroscopic properties. 

Measurement of the SN la rate may also ha ve a more direct i mpac t on the determination 
of syst e matic error s in SN la d istance estimates. iMannucci et al.l (120061 ) ; IScannapieco fc Bildsten 
(2005); Neill et al. (2006) and Sullivan et al. ( 20061 ) have argued that a two-component model 
of the SN la rate, in which a prompt SN component follows the star formation rate and a 
second component follows the total stellar mass, is strongly favored over a single SN la 
channel. In this picture, since the cosmological star formation rate increases sharply with 
lookback time, the promp t comp onent is expected to dona inate the total SN la rate at high 
redshift. IMannucci et al.l (120061 ) and iHowell et al.l (120071 ) pointed out that this evolution 
with redshift can be a potential source of systematic error in SN la distance estimates, if the 
two populations have different properties. 

In order to test such a model for the evolution of the SN la rate, improved measurements 
of the rate as a function of redshift and of host galaxy properties are needed. The Supernova 
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Legacy Survey (SNLS) has recently presented th e most precise naeasurement of the S N la 



rate at high redshift {z ~ 0.5) based on 73 SNe la fiNeill et al.ll2006 
low redshifts (z ~ 0.1) , SN la rate measu rements (ICappellaro et al. 



Sullivan et al. 12006). At 



1999 



Hardin et al. 



2000 



Madgwick et al.ll2003l : iBlanc et al.ll2004j ) have suffered fron a small sample sizes and also from 
systematic errors associated with heterogeneous samples (ICappellaro et al.lll999l) and with 



selection biases due to the targeting of known, relatively l uminous galaxies ( Hardin et al. 
2000l : iBlanc et al.l 12004 ) . The low-redshift measurement of iMadgwick et al.l (120031 ) , based 
on SNe discovered fortuitously in SDSS galaxy spectra, is affected by different systematic 
uncertainties than traditional photometric searches, e.g., due to the finite aperture of the 
SDSS spectroscopic fibers. 

In this paper, we present a new measurement of the SN la rate at low redshift, based 
upon th e first season of data from the SDSS-II Supernova Survey. The SDSS-II Supernova 
Survey (IFrieman et al.l l2008l ) offers several advantages for this measurement. It covers a 
larger spatial volume than previous SN surveys, a result of the combination of intermediate- 
scale (2.5-m) telescope aperture, wide field of view (3 square degrees), modest effective 
sidereal exposure time (54 sec), and use of drift- scanning to efficiently cover a large sky 
area (~ 300 square degrees). The SDSS-II Supernova Survey is a rolling search, with new 
SNe discovered simultaneously with the follow-up of previously discovered SNe. Unlike SN 
searches that target known galaxies, the SDSS-II Supernova Survey is not biased against 
finding SNe in low-luminosity ho st galaxies. Well-calibrated photometry in the SDSS ugriz 
passbands (IFukugita et al.lll996l ). with a typical interval between observations of four days, 
yields well-sampled, multi-band light curves that enable photometric typing of SNe with 
high confidence. Moreover, rapid on-mountain photometric reduction and image processing 
coupled with an extensive spectroscopic follow-up program enable spectroscopic confirmation 
of a very high fraction of the low-redshift SN la candidates. 

The SDSS-II Supernova Survey was carried out over three three-month seasons, during 
Sept. -Nov. 2005-7. The results presented here are based on the Fall 2005 season. The SDSS- 
II Supernova Survey measures light curves for SNe la to redshift z ~ 0.4, with a median 
redshift of (z) = 0.22 for spectroscopically confirmed SNe la. However, in this first paper we 
limit the analysis to low redshift, z < 0.12, since our spectroscopic follow-up is essentially 
complete over this redshift range, and the uncertainty due to spectroscopically unobserved 
(and untyped) SNe is therefore negligible. In presenting a SN la rate measurement, one must 
decide whether to include pecuhar SNe la, i.e., events that are photometrically and/or spec- 
troscopically unusual, since it is not clear that they are members of the same population as 
the "normal" SNe la. Formerly, the peculiar designation included events such as 1991T and 
1991bg, which are highly overluminous and underluminous events, respectively. However, 
since these SNe appear to follow the standard peak-luminosity/decline-rate relation, they are 
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now generally considered extreme n iembers of the norm al SN la po pulation (INugent et al. 
1995[ ). Other events, such as 2002ic jHamuv et al-lboosh and 2002cx Jli et alj[2003l ). exhibit 
more pronounced peculiarities and do not fit the luminosity-decline relation. The first season 
of the S pSS-II Supernova S urvey includ ed two such truly peculiar events at lo w redshift, 
2005hk JphiUips et aliboOTh and 2005gj JPrieto et aliboOTl : kldering et al.lbooeh . Although 
these peculiar events may arise from the same evolutionary path as normal SNe la, which 
would argue for including them in a SN la rate measurement, we have chosen to include only 
SNe with light curves that obey the standard brightness-decline relation. More specifically, 
we include in our rate measurement sample only SNe with light curves that are well described 
by the MLCS2k2 SN la light curve model jRiess et al.lll99d : Ijha et all 120071 ) . see Re- 
gardless of the physical arguments surrounding peculiar events, we exclude them from this 
analysis primarily because we do not yet have a robust determination of our efficiency for 
detecting them. 

The rest of the paper is organized as follows. In §2] we provide a brief outline of the 
survey observing strategy and operations as they relate to the rate determination. In ^we 
define selection criteria and present the sample of SNe la used in this measurement, based 
on spectroscopic and photometric measurements. In ^we present estimates of the detection 
efficiency for low-redshift SNe la, based on artificial SNe inserted into the survey images and 
on Monte Carlo simulations. We present our measurement of the SN la rate and discuss the 
SN la rate as a function of host galaxy type in §5] In ^ we compare our result to other SN 
la rate measurements and combine rate measurements to fit semi-empirical models of rate 
evolution. 



2. SDSS-II Supernova Survey Overview 



Here we briefiy describe the SDSS-II Supernova Survey, highlighting the fe atures that are 



d 



m 



most relevant to a rate measurem ent. The survey is described in more detail in lFrieman et al. 
f lgoggh and 



Sako et al.l (120081 ). A tech nical summary of t he SPSS is given b y 



20001) . and further details c an be found iiilHogg et al.l (120011 ): llvezic et al.l (12004 ) : lLupton et al 



1999[ ): ISmith et all J20021 ): iTucker et al.N2006h . 



York et al. 



2.1. Imaging 



The SDSS-II Supernova Survey is carried out on the 2.5m telescope (|Gunn et al. 



at Apache Point Observatory (APO), using a wide- field CCD camera (IGunn et al. 



(Gunn et al. 


2006) 


(IGunn et al. 


19981) 
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operating in time-delay-and-integrate (TDI, or drift sc an) mode. Observati ons are obtained 



nearly simultaneously in the SDSS ugriz filter bands (IFukugita et al.lll996l ) 



The SDSS-II Supernova Survey covers a region, designated stripe 82, centered on the 
celestial equator in the Southern Galactic hemisphere, bounded by —60° < aj2ooo < 60°, 
and —1.258° < Sj2ooo < 1.258°. Stripe 82 has been imaged multiple times in photometric 
conditions by the SDSS-I survey; co-added images from those runs provide deep template 
images and veto catalogs of variable objects for the SDSS-II Supernova Survey transient 
search. Due to gaps between the CCD columns on the camera, each stripe is divided into 
northern (N) and southern (S) strips; the SDSS-II Supernova Survey alternates between the 
N and S strips on subsequent nights. Each strip encompasses ~ 162 square degrees of sky, 
with a small overlap between them, so that the survey covers ~ 300 square degrees. On 
average each part of the survey region was observed once every four nights during the 2005 
season. Figure [1] shows the sky coverage versus survey time, along with a representative 
SN la light curve. 
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Fig. 1. — Right ascension range covered by SDSS-II Supernova Survey imaging runs 
vs. epoch. The panels labelled N & S denote the Northern and Southern strips of stripe 
82. The regions aj2ooo < —51° and aj2ooo > 57° are not covered early in the season, and 
these regions are suppressed from the rate measurement. The top panel shows an example 
unextincted g-band light curve for a SN la at a redshift of 0.12, based on the MLCS2k2 
model. 
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2.2. Supernova Search Pipeline 

There are five main components to the supernova search pipeline: photometric reduc- 
tion, image subtraction, automated object selection, visual inspection, and light curve fitting 
for spectroscopic target selection. We describe them briefiy in turn. For a full night of imag- 
ing data, the entire pipeline runs in approximately 20 hours, sufficient for keeping up with 
the data fiow and for rapid spectroscopic targeting. 

In the first stage of the search pipeline, the imaging data is acquired from the camera and 



proce ssed through the the SDSS photometric reduction pipeline, known as PHOTO (iLupton et al 



2001 ). PHOTO produces "corrected" images that are astrometrically calibrated (iPier et al. 
20031 ) and provides a local estimate of the point spread function (PSF). In the second stage, 
images are processed through the difference imaging pipeline. To run the search pipeline 
to completion in less than a day with the available on-mountain computing resources, only 
the corrected gri images are processed beyond the first stage. The search image is astro- 
metrically and photometrically registered to the template image, an d the template image i s 



convolved with a kernel chosen to minimize subtraction residuals (lAlard fc Luptod Il998l ). 
A difference image is then obtained by subtracting the convolved template image from the 
survey image. Peaks are d etected in the differen ce image using the DoPHOT photometry and 



object detection package (ISchechter et al.lll993l ). The signal-to-noise threshold for object 



detection is at ~ 3.5, corresponding in typical conditions to (7 ~ 23.2, r ~ 22.8, and i ~ 22.5. 
The typical magnitudes at signal-to-noise of 10 for point-like objects are g ~ 21.8, r ~ 21.5, 
and i ~ 21.2. 

The third stage of the SN search pipeline comprises a sequence of automated filtering 
operations that select events of potential interest from among those detected in the differ- 
ence images. We require a close positional match in at least two of the gri images, which 
removes cosmic rays, single-band spurious noise fluctuations, and a large fraction of asteroids 
and other rapidly moving objects detected by the survey. All detections that satisfy these 
criteria are entered into a MySQL database and are referred to as objects. To reject active 
galactic nuclei (AGN) and variable stars, we veto any detection occurring at the position of 
a previously cataloged variable, using observations of stripe 82 from several previous years. 
The area corresponding to previously cataloged variable objects represents ~ 1% of the total 
survey area. 

In addition to SNe, the database of detected objects includes a variety of physical and 
non-physical transients. Physical sources include slow-moving asteroids that were not re- 
jected by the moving object veto, AGN and variable stars not already cataloged, and high 
proper-motion stars. Non-physical sources include improperly masked diffraction spikes from 
bright stars and artifacts of imperfect image registration. To remove non-physical sources. 
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cut-out images of all objects that remain after the automated filtering are visually inspected 
and classified in the fourth stage of the search pipeline. Objects visually classified as con- 
sistent with a possible SN event are flagged for further analysis and are denoted candidates. 
Subsequent object detections in difference images at the same position are automatically 
associated with the same candidate. 

In the fifth and final stage of processing for the SN search, the gri light curve for each SN 
candidate is fit to models of type-la, type-Ib/c and type-II SNe. The non-la SN models con- 
sist of template light curves constructed from photometric measurements of individual SNe 
provided b y the SUSPECT dat abasjll, coupled with the corresponding SN spectral model 



provided by lNugent et al.l (120021 ). For the SN la model, a stretch and a wavelength-dependent 



scale factor is applied to a fiducial bolometric light curve in a way de signed to reproduc e 



the Ami5 parameterization of the peak-luminosity /decline-rate relation (iHamuy et al.lll996l ). 
The time of maximum. Amis, redshift, and extinction parameter Ay (magnitudes of extinc- 
tion in the V-band) are fit parameters that are searched on a grid for the set of values that 
produce the minimum value of the statistic. For some candidates, we additionally carry 
out difference imaging in the u and z passbands in order to better distinguish Type II and 
Type la SNe that tend to have a significantly different u — g color at early epochs. To further 
constrain the early light curve shape, we carry out forced-positional photometry on difference 
images at the position of the candidate in pre-discovery images. The relative goodness of fit 
of candidate gri hght curves to SNe la and core-collapse SNe models is used as a factor in 
prioritizing spectroscopic follow-up. In particular, all SN la candidates found before peak 
and with estimated current r-band magnitude < 20 are placed on the spectroscopic target 
list, and our follow-up observations are nearly complete out to that magnitude. Since the 
typical peak magnitude for a SN la with no extinction at redshift z = 0.1 is r ~ 19.3, we 
might expect that the spectroscopic SN la sample should be essentially complete out to 
roughly this redshift as well; we shall see later that this is the case. This photometric pre- 
selection of SNe la has proven very effective: approximately 90% of the candidates initially 
targeted as SNe la after two or more epochs of imaging have resulted in a SN la spectroscopic 
confirmation. The SDSS-II Supernova Surve y photometric cla ssification and spectroscopic 



target selection are discussed in full detail in ISako et al.l (120081 ). 



^http://bruford. nhn.ou.edu/suspect/indexl. html 
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2.3. Artificial Supernovae 

To measure the SN rate, it is clearly important to understand the efficiency of the survey 
for discovering SNe. As part of normal survey operations, we insert artificial SNe la (hereafter 
fakes) directly into the corrected survey images after the photometric reduction (PHOTO) but 
before difference imaging. The primary motivation for inserting fakes into the data stream is 
to provide real-time monitoring of the performance of the survey software pipeline and of the 
human scanning of objects. The fakes provide quantitative information about the efficiency 
of the survey software, human scanning, and the photometric classification of SNe la that is 
useful in the rate determination. Here we describe the basic algorithm for gene rating fakes 



and inserting them into the data stream; for more details, see ISako et al.l (120081 ). 



A fake is a pixel-level simulation of a point source with a light curve chosen to closely 
represent that of a real SN la. At each epoch for which the fake has a chance of being 
detected, the calculated CCD signal for the fake is directly added to the survey image. For 
the 2005 observing season, we generated a library of 874 fake light curves: each fake light 
curve is assigned a position, redshift, date of peak luminosity (in V-band), and an intrinsic 
luminosity that correlates with decline rate. This resulted in ~ 7, 800 fake epochs during 
the season. The redshift distribution for the fakes was generated by assuming that the 
number of SNe la is roughly proportional to the volume element, {dN/dz) oc z^, in the range 
0.0 < ;z < 0.4. 

To model the effect of contamination from host galaxy light on the detect ion efficiency, 



each fake is placed near a galaxy selected from the photometric redshift catalog (lOyaizu et al. 



20071 ) for SDSS imaging on stripe 82. A host galaxy is drawn at random, from a distribution 
proportional to the r-band luminosity, from galaxies which have a photometric redshift within 
~ 0.01 of the redshift assigned to the fake. 

The SN la light curve model used to generate ugriz magnitudes for a fake at each epoch 
is the same model that is used for early light curve ffiting and photometric typing on the 
imaging data, but with the light curve parameters now chosen from an input probability 
distribution. To generate a point-source image from the ideal magnitudes, we use the es- 
timate of the PSF from PHOTO at the position of the fake at the given epoch. We obtain 
the conversion from magnitudes to instrumental units (analog-to-digital units, or ADU) by 
running the DoPHOT photometry package on a set of cataloged stars in the survey image for 
which the magnitudes have been previously measured by the SDSS. After scaling the PSF 
model to match the computed ADU flux, we add Poisson fluctuations to each pixel. Finally, 
the row and column in the field that correspond to the position of the fake are taken from 
the astrometric solution provided by the imaging pipeline, and the fake is overlaid on the 
survey image. 
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When a fake is detected in the difference images, its identity as a fake is kept hidden 
while it is scanned by humans. After scanning, the fakes are revealed so that they are 
not mistakenly targeted for spectroscopic follow-up and so that the efficiency of scanners in 
tagging fakes as SN candidates can be monitored. However, like all candidates, fakes are 
processed through the automated hght curve fitter /photometric typing algorithm so that we 
can test if they are accurately typed as SNe la after a few photometric epochs. The use of 
the fakes for measuring the survey detection efficiency is discussed in § 14.11 



2.4. Spectroscopy 



The classification of SNe is defined by their spectroscopic features. In addition, spec- 
troscopy provides a precise redshift determination and, in a number of cases, host galaxy 
spectroscopic-type information. Spectroscopic follow-up of the SDSS-II Supernova Sur- 
vey candidates is being undertaken by a number of telescopes. During the 2005 observing 
season, spectroscopic observations were provided by the Hobby-Eberly 9.2m at McDonald 
Observatory, the Astrophysical Research Consortium 3.5m at Apache Point Observatory, the 
William-Herschel 4.2m, the Hiltner 2.4m at the MDM Observatory, the Subaru 8.2m and 
Keck 10m on Mauna Kea, and the SALT 11m at the South African Astronomical Observa- 
tory. 

The classification of SN spectra is performed by compari ng the spectra l data to normal 
and peculiar supernova spectral templates from th e work of iNueent et al.l (120021) and to a 



publi c library of well-measured supernova spectra (IMatheson et al 



12005 



Blondin fc Tonry 



20071 ). The SN typing in this work is based on vis ual inspection of the spectra, but was 
guided by applying the cross-correlation technique of iTonry fc David (Il979l ) to the spectrum 
and the template library. The visual inspection relies heavily on the characteristic SN la 
features of Si and S absorption, which are usually prominent at optical wavelengths for this 
redshift range. 

The redshift determination is based on galaxy features when they are present; otherwise 
SN features are used. In some cases, particularly at low redshift, a high-quality spectrum of 
the SN host galaxy is available from the SDSS-I spectroscopic survey. Comparison with those 
spectra indicate that our follow-up spectroscopic redshifts are determined to an accuracy of 
~ 0.0005 when galaxy features are used and ~ 0.005 when SN features are used. Further 
de tails of the SDSS- II Supernova Survey first-season spectroscopic analysis are presented 
in 



Zheng et all fl2008f ) 
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2.5. Final Photometry 

To obtain more precise SN photometry than the on-mountain difference imaging pipeline 
provides, we re-process the imaging data for all spectroscopi cally confirmed and ot her inter- 



esting SN candidates through a final photometry pipeline (iHoltzman et al.ll2008l ). In this 
"scene-modeling photometry" (SMP) pipeline, the supernova and the host galaxy (the scene) 
are modeled respectively as a time-varying point-source and a background that is constant 
in time, both convolved with a time- varying PSF. This model is constrained by jointly fitting 
all available images at the SN position, including images well before and after the SN explo- 
sion. Since there is no spatial resampling or convolution of the images that would correlate 
neighboring pixels, the error on the flux can be robustly determined. The SMP pipeline often 
provides photometric measurements at additional epochs compared to the survey operations 
pipeline. The final analysis of SN light curves discussed in this paper is based on SMP; in 
particular, the selection cuts described in ^ are made using the SMP pipeline. 



3. Defining the SN la Sample for the Rate Measurement 

The SN la sample for the rate measurement must include all SNe la in the redshift 
range of interest, not just those for which we have a confirming spectrum. Although we have 
high efficiency for discovering and spectroscopically confirming low- redshift SNe la ( §3.ip . 
we can take advantage of our rolling search data to carry out an extensive post-season hunt 
for SNe la that may have been missed by the search pipeline during the survey season ( §3.2p . 



3.1. Spectroscopic SN Sample 

In its first season (Fall 2005), the SDSS-II Supernova Survey discovered 130 events with 
secure spectroscopic identifications as SNe Icfl and 16 events that are considered probable 
SNe la based on their spectra. For SN la events satisfying the selection criteria below, the 
spectroscopic follow-up is essentially complete for redshifts z < 0.12, so we have chosen 
to focus on this redshift range for a first measurement of the SN rate. For z < 0.12, the 
sample contains 27 spectroscopically confirmed SNe la and 2 spectroscopically probable SNe 
la before making selection cuts. 

For the measurement of the SN la rate, we impose a number of selection criteria on 



^The classification of 2005gj as a SN la may be controversial (jPrieto et al. 
exclude it from this analysis. 



20071) : as noted in gH we 
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the SN photometric data, with the aim of producing a sample that has a well-characterized 
selection function. These criteria are applied to the spectroscopically confirmed and prob- 
able SNe la with z < 0.12. For consistency, we will also apply these selection cuts to the 
photometric (i.e., spectroscopically unconfirmed) SN sample discussed in §3.2[ The selection 
cuts for the rate measurement are as follows: 



1. -51° < «j2ooo < 57°. 

Although the SDSS-II Supernova Survey covers the RA range —60° < aj2ooo < 60°, 
early in the Fall 2005 observing season we did not have complete templates available 
for the regions aj2ooo < —51° and aj2ooo > 57°, so these RA regions were not initially 
used for the SN search, as shown in Fig. [H In principle, we could account for this 
by modeling the time-varying effective search area, but for simplicity we choose to 
excise these RA regions from the rate measur ement. Furthermo re, the calibration star 



catalog used by our final photometry pipehne (llvezic et al.ll2007l ) does not extend below 
'^j2ooo ~ "~51°, and we cannot presently simulate light curves for SNe in this region. 
This cut removes one confirmed SN la, 2005iu, from the rate sample. 

2. There are photometric observations on at least five separate epochs between —20 days 
and +60 days relative to peak light in the SN rest-frame. 

Peak light refers to the date of maximum luminosity in the SN rest-frame B-band 
according to the best-fit MLCS2k2 light curve model. This cut requires that the 
light curve is reasonably well-sampled and it is primarily useful for photometrically 
distinguishing SNe la from other SN types with high confidence when there is no SN 
spectrum available (see §3.21) . Here and below, a photometric observation simply means 
that the survey took imaging data at that epoch on that region of sky and that SMP 
reported a SN flux measuremen t (not necessarily significant or even positive) with no 



error fiags (see iHoltzman et al.l (120081 )) in at least one of the three gri passbands. It 
does not imply a detection above some signal-to-noise threshold. One SN discovered 
late in the observing season, 20051k, fails this cut. 

3. At least one epoch with signal-to-noise ratio > 5 in each of g, r, and i (not necessarily 
the same epoch in each passband). 

This cut ensures that there are well-measured points on the light curve, and is mainly 
useful for rejecting low signal-to-noise events from the photometric sample. All spec- 
troscopically confirmed SNe la in the low-redshift sample satisfy this cut. 



4. At least one photometric observation at least two days before peak light in the SN rest 
frame. 
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5. At least one photometric observation at least ten days after peak light in the SN rest 
frame. 

These two cuts require sampling of the light curve before and after peak light, ensuring 
that we have a precise determination of the time of peak light. These cuts also help 
remove non-SN la contaminants from the photometric sample (see §3.2p . Finally, they 
guarantee that the epoch of peak light occurs during our observing season, i.e., between 
Sept. 1 and Nov. 30, which is one of the criteria used in defining the rate measurement 
in ^ Since these cuts are more restrictive than the requirement of peak light during 
the observing season, they are the main contributors to the inefficiency estimated in 
§4.21 These are the most restrictive cuts on the spectroscopic SN la sample, together 
removing nine events: four SNe la discovered early in the observing season do not have 
a pre-maximum observation, and five SNe la found late in the season do not have a 
photometric observation more than ten days past peak light. 

6. MLCS2k2 light curve fit probability > 0.01. 



The MLCS2k2 light curve fitter flRiess et al.lll996l : Ijha et al.l 120071 ) takes as input the 



measured SN magnitudes in each passband at each epoch, and the measured SN red- 
shift; it then finds the likelihood as a function of the four parameters /i (the distance 
modulus), Av (the extinction parameter), the time of peak light in rest-frame B-band, 
and the light curve shape/luminosity parameter A. The MLCS2k2 fit probability is 
defined by evaluating the usual statistic for the data and the best fitting MLCS 
model and assuming that this statistic obeys a Xn-4 probability distribution, where n 
is the number of photometric data points. The model parameters of the best fitting 
MLCS2k2 model are defined as the mean of the probability distribution for each corre- 
sponding parameter. This cut on the fit probability provides an automated method of 
removing photometrically peculiar SNe la from the sample. We find that essentially all 
of the spectroscopically normal SNe la in our confirmed sample have a fit probability 
> 0.1. However, the spectroscopically confirmed SN la sample is likely to be biased 
toward "high-quality" light curves, so we place the selection cut at a less restrictive 
value. Three spectroscopically identified SNe la are rejected by this cut, including 
the peculiar SNe 2005hk (xVd.f. = 90/21) and 2005gj (xV^-f- = 198/45). The third 
rejected SN, with internal SDSS candidate designation 6968 (x^/d.f. = 78/27), was 
classified as a spectroscopically probable SN la (see §2.41) and shows some evidence of 
being spectroscopically similar to 2005hk. For the sample of photometric SN candi- 
dates ( §3.2p . this cut helps remove non-SN astrophysical variables, such as AGN and 
M stars. 



7. MLCS2k2 light curve fit parameter A > -0.4. 



- 16 - 



The MLCS parameter A is a measure of the hght curve shape and intrinsic luminosity. 
Smaller values of the A parameter correspond to more slowly-declining, intrinsically 
brighter SNe la. This cut requires that A be consistent with the values observed for 
the low-redshift SNe la that were used to train the MLCS2k2 light curve fitter. For 
the photometric SN candidates, this cut helps reject Type II supernovae, which often 
have a long plateau after the epoch of peak luminosity and result in large negative 
fitted values of A. 

The above selection requirements result in the 16 spectroscopically identified SNe la 
that are listed in Table [H This sample includes 2005je, which was classified as spectro- 
scopically probable. The spectroscopically confirmed SNe la that are removed from the 
rate-measurement sample are listed in Table [21 the last column indicates which of the above 
selection criteria was used to remove each SN. 



Table 1. SNe la included in the rate sample. 



SDSS Id lAUC a (J2000.0) 5 (J2000.0) Redshift Redshift 
Designation Source 



1 241 


2005fF 


22 


30 


41 


41 


—00 

WW 


46 

t:W 


35 


7 


n 

W 


088 

WOO 


SN 


1371 


2005fh 


23 


17 


29 


71 


+00 


25 


45 


8 





120 


galaxy 


2561 


2005fv 


03 


05 


22 


42 


+00 


51 


30 


1 





119 


galaxy 


3256 


20051in 


21 


57 


04 


23 


-00 


13 


24 


4 





107 


galaxy 


3592 


2005gb 


01 


16 


12 


58 


+00 


47 


31 








086 


galaxy 


3901 


2005ho 


00 


59 


24 


10 


+00 


00 


09 


3 





063 


galaxy 


5395 


2005hr 


03 


18 


33 


81 


+00 


07 


24 


3 





117 


SN 


5549 


2005hx 


00 


13 


00 


13 


+00 


14 


53 


7 





120 


SN 


5944 


2005hc 


01 


56 


47 


94 


-00 


12 


49 


1 





046 


galaxy 


6057 


2005if 


03 


30 


12 


87 


-00 


58 


28 


5 





067 


galaxy 


6295 


2005js 


01 


34 


41 


51 


-00 


36 


19 


4 





084 


SN 


6558 


20051ij^ 


01 


26 


48 


40 


-01 


14 


17 


3 





057 




6962 


2005je 


02 


35 


26 


61 


+01 


04 


29 


6 





094 


galaxy 


7147 


2005jli 


23 


20 


04 


42 


-00 


03 


19 


8 





109 


galaxy 


7876 


2005ir 


01 


16 


43 


80 


+00 


47 


40 


7 





076 


galaxy 


8719 


2005kp 


00 


30 


53 


15 


-00 


43 


07 


9 





117 


galaxy 


9266^ 




03 


20 


43 


16 


-01 


00 


07 


2 





036 


galaxy 



^SN type confirmed by lQuimby et al.l (120051 ) 



^'Photometrically identified SN la. See § (l3.2p 

Note. — SDSS Id denotes internal candidate designation. 
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Table 2. Spectroscopically confirmed SNe la with z < 0.12 cut from the rate sample. 



SDSS Id lAUC a (J2000.0) 6 (J2000.0) Redshift Cut Index 
Designation 



722 


2005ed 


00 


02 


49.37 


+00 


45 


04.6 


0.086 


4 


739 


2005ef 


00 


58 


22.87 


+00 


40 


44.6 


0.107 


4 


774 


2005ex 


01 


41 


51.24 


-00 


52 


35.0 


0.093 


4 


2102 


2005fn 


20 


48 


53.04 


+00 


11 


28.1 


0.095 


4 


4524 


2005gj 


03 


01 


11.95 


-00 


33 


13.9 


0.062 


6 


6773 


2005iu 


20 


20 


15.61 


+00 


13 


02.5 


0.090 


1 


6968 




01 


18 


13.37 


-00 


54 


23.6 


0.098 


6 


8151 


2005hk 


00 


27 


50.88 


-01 


11 


53.3 


0.013 


6 


10028 


2005kt 


01 


10 


58.04 


+00 


16 


34.1 


0.066 


5 


10096 


20051j 


01 


57 


43.03 


-00 


10 


46.0 


0.078 


5 


10434 


20051k 


21 


59 


49.43 


-01 


11 


37.3 


0.103 


5,2 


10805 


2005ku 


22 


59 


42.61 


-00 


00 


49.3 


0.045 


5 


11067 


2005ml 


02 


14 


04.42 


-00 


14 


21.1 


0.119 


5 



Note. — SDSS Id denotes internal candidate designation. See section 13.11 
for explanation of cut index. 
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3.2. Photometric SN Sample 

In addition to tlie spectroscopically identified SNe la discussed above, the survey has 
measured light-curves for a few thousand variable objects, including possible SNe, for which 
we did not obtain a classifiable spectrum while the source was bright enough to identify. 
We refer to these spectroscopically unobserved or unclassified objects as photometric SN 
candidates. There are a number of reasons for such spectroscopic incompleteness, including 
limited spectroscopic resources, targeting errors (e.g., misplacement of a spectroscopic slit), 
poor weather either preventing spectroscopic observations or rendering them indeterminate, 
and possible inefficiencies in the spectroscopic target selection algorithm. In order to make a 
reliable SN la rate measurement, we must investigate the photometric SN candidates to 
determine the level of incompleteness, if any, of the spectroscopic SN la sample. This 
is a challenge, because a sample of purely photometric SN la candidates may be heavily 
contaminated by objects that are not SNe la, especially if there are significant numbers of 
objects with multi-band light curves that are not too dissimilar from those of SNe la. The 
combination of selection cuts listed in §3.11 is designed to meet this challenge, by rejecting 
the majority of non-SN la contaminants. In addition to the spectroscopically confirmed and 
spectroscopically probable SNe la discussed in §3.1, the SDSS-II Supernova Survey discovered 
16 low-redshift SNe that were spectroscopically confirmed as non-la SNe in it's first year. 
As a check that our selection cuts are effective at rejecting non-la SNe, we apply the same 
cuts to this sample of 16 low-redshift {z < 0.2) spectroscopically confirmed non-la SNe. All 
but one of these non-la SNe are rejected by these selection cuts. The selection criteria above 
could be made more restrictive in order to reduce potential non-la contamination of the 
photometric SN sample. For example, by requiring a photometric observation at least 16 (as 
opposed to 10) days after peak light in the SN rest frame, the spectroscopic non-la SN above 
would be eliminated from the sample. However, we find that such a change would have no 
impact on the selection of photometric SN candidates for inclusion in the rate sample. 

To determine whether any of the photometric SN candidates are genuine SNe la in the 
redshift range of this rate measurement, we must estimate both the SN type and redshift for 
each candidate. There are two categories of photometric SN candidates, (a) those for which 
we have a precise spectroscopic measurement of the redshift and (b) those for which we do 
not. The redshifts for category (a) candidates come from two sources. The first source is 
the SDSS-I spectroscopic galaxy survey, which measured redshifts for ~ 28, 000 galaxies in 
our survey region at redshifts z < 0.12. The second source is from subsequent spectroscopic 
observations of ~ 80 host galaxies of the highest-quality photometric SN candidates; these 
spectra were obtained in the summer and fall of 2006 and 2007. Using the sample selection 
process described below in §3.2.21 we found in our imaging data only one photometric SN la 
candidate that passes the selection criteria in §3.11 and that has a spectroscopic redshift 
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z < 0.12 (category (a)). The host galaxy of this SN la candidate, which has internal SDSS 
SN designation 9266, has a spectroscopic redshift of z = 0.0361 measured by the SDSS galaxy 
redshift survey. This object was not targeted for spectroscopic follow-up during the SDSS-II 
Supernova Survey because it has very high extinction. Ay ~ 4 according to the MLCS2k2 
fit. This extinction value lies outside the range of the Ay-search grid for the photometric 
typing algorithm used during the search ( §2.21) . 



3.2.1. Redshift estimation for photometric SN candidates 

For each photometric SN candidate without a spectroscopically determined redshift (cat- 
egory (b) above), we must estimate both the redshift and the SN type from the photometric 
data. We do this using a modification of the standard MLCS2k2 light curve fit, in which the 
redshift is included as a parameter in the likelihood function. In this instance, the distance 
modulus /i is not treated as a fit parameter; instead, we adopt the concordance LCDM cos- 
mology, with VLjn = 0.3, = 0.7, and dark energy equation of state parameter w = — 1, 
and fix fi{z) to its functional form for that cosmology. The photometric redshift estimate 
Zphot is then obtained by marginalizing over the other fit parameters, i.e., the epoch of peak 
luminosity, the extinction Ay, and the shape parameter A. 

Although these SN photometric redshift estimates depend on the assumed cosmology, 
we do not expect them to be extremely sensitive to the values of the cosmological param- 
eters, especially at the modest redshifts under consideration here. To test the accuracy of 
these redshift estimates, we applied the MLCS2k2 redshift fit to SNe la light curves that 
have spectroscopically measured redshifts. For this test, we use two sets of objects: (i) all 
spectroscopically confirmed SNe la with redshift z < 0.25, and (ii) all photometric SN la 
candidates that satisfy the selection criteria in §3. II and that have spectroscopic (host galaxy) 
redshift z < 0.25. We include objects in category (ii) because the spectroscopically confirmed 
SNe la could represent a biased sample of the SN la population if there are spectroscopic 
selection biases. We include objects with redshift z > 0.12 to yield a more statistically 
significant test and to check for photometric redshift biases that could cause these objects 
to be erroneously included in the z < 0.12 sample. There are 61 and 28 events in categories 
(i) and (ii), respectively. 
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Fig. 2. — Residuals of the photometric redshift estimates, Zphot 
sample of spectroscopically confirmed SNe la (black points) and for the photometric SN la 
candidates that satisfy the rate selection cuts and for which host galaxy redshifts are available 
(open points). The marginahzed redshift errors reported by the MLCS2k2 hght-curve fits 
are shown. 
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Fig. 3. — The distribution of photometric redshift residuals for the spectroscopically con- 
firmed SNe la (left panels) and the photometric SN la candidates (right panels) shown in Fig. 
El Upper panels show distributions of the difference between the photometric redshift, Zphot, 
and the spectroscopic redshift, Zspec] lower panels show distributions of {zphot — Zspec)/o'z, 
where is the photometric redshift uncertainty reported by the MLCS2k2 fit. Inset panels 
show the inferred mean and dispersion of the Gaussian fits to each distribution. 
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Fig. 4. — Left panel: The distribution of photometric redshift residuals, Zphot — Zgpec, for 
the combined spectroscopic+photomctric SN la samples. Zp^ot is the photometric redshift 
and Zgpec the spectroscopic redshift; Right panel: distribution of residuals normalized by the 
reported photometric redshift uncertainty, a^- 
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The residuals of the SN photometric redshift estimates for this test sample of 89 objects 
are shown in Figure [2l The distributions of the residuals are shown in Figure [3] separately 
for categories (i) and (ii), which indicates that the distributions for the two samples are 
consistent. The fit residuals for the combined sample are shown in Figure |H For the combined 
sample, the mean residual of the photometric redshift estimate is consistent with zero. The 
scatter in the photometric redshift estimate for the combined test sample is az = 0.018, 
and Fig. [2] shows that the scatter increases with redshift. The distribution of the residuals 
normalized by the MLCS-reported redshift error is shown in the right panel of Figure |H If the 
reported redshift errors were accurate and Gaussian, this distribution would be a Gaussian 
with unit variance, cr = 1. The distribution appears to be approximately Gaussian, but with 
measured variance a = 1.4; we therefore choose to multiply the MLCS-estimated photometric 
redshift error for each candidate by 1.4. 

In addition to the SN photometric redshift estimates, we also have hos t galaxy pho 



tometric redshift estimates for the majority of photometric SN candidates (jOyaizu et al. 



20071 ). Although the galaxy Zphot estimates have larger scatter than those from the SN light 
curves, in principle we could require consistency between these two redshift estimates as an 
additional selection cut on the photometric SN sample. Since core-collapse SNe are typically 
fainter than SNe la, they would typically be assigned incorrectly high photometric redshifts 
by the light curve fitter. Using the existing selection cuts, however, we find no contamination 
of the rate-measurement sample from the photometric SN candidates without spectroscopic 
redshifts (see §3.2.21) . Therefore a requirement of consistency between the supernova- and 
galaxy-derived photometric redshift estimates is not necessary in the present analysis. 



3.2.2. Selection of photometric SN candidates 

In the Fall 2005 observing season, the software and human data processing pipeline 
described in §2.2l yielded 11,385 SN candidates, including the 146 spectroscopically confirmed 
and probable SNe la and 20 that were confirmed as other SN types. The majority of the 
remaining candidates (~ 60%) are single epoch events that are most likely to be slow- 
moving asteroids, leaving ~ 4500 multi-epoch SN candidates. To search for photometric 
SNe la among this large set of candidates, we studied two subsamples selected according to 
different criteria. 

The first photometric subsample is designed to exhaust the list of candidates that are 
most likely to be SNe la. This subsample was selected by choosing all SN candidates that 
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the survey photometric typing code (described in §2.21) classified as SNe Ic@ and that were 
detected in at least 3 epochs by the on-mountain software pipeline. The images for the 
resulting subsample of ~ 420 candidates were processed through the final SMP pipeline, 
and the resulting light curves were fit with the MLCS2k2 program, using the redshift as a 
fit parameter in cases where there was no measured host galaxy spectroscopic redshift (the 
majority of cases). One highly extincted SN la (SDSS-SN 9266, discussed in §3.21) . with a 
host galaxy redshift measured by the SDSS galaxy survey, was recovered from this subset of 
photometric SN candidates. No other candidates in this subsample pass the rate selection 
cuts and have a spectroscopic or SN photometric redshift z < 0.12. 

The second photometric subsample is designed to study the candidates that are less 
likely to be SNe la. This subsample was selected by choosing all SN candidates with detec- 
tions at more than two epochs during the search and with an estimated time of maximum 
light, based on the survey photometric typing code, in the twenty-day interval between mod- 
ified Julian day (MJD) of 53660 and 53680 (October 17 and November 6). Since the selection 
criteria for this second subsample are looser than for the first (there is no requirement that a 
SN la light-curve template provides the best-fit), the number of candidates it selects would 
be an order of magnitude larger. Using a restricted time interval provides a manageable 
number of events to study that are representative of the population of these lower-quality 
light curves. These selection criteria result in 462 candidates, which represent ~ 1/6 of the 
multi-epoch SN candidates that have not already been included in the samples discussed 
above. We find no events in the second subsample that pass the rate measurement selection 
cuts and that have a spectroscopic or SN photometric redshift z < 0.12. 

Although no other photometric SN candidates pass our selection criteria, we must al- 
low for uncertainties in the SN photometric redshift estimates from MLCS2k2. In the two 
photometric subsamples, two candidates that pass the rate-selection cuts have estimated SN 
photometric redshifts within ~ 1.5cr of our cutoff of z = 0.12, using the inflated redshift 
errors discussed in §3.2.11 One of these candidates is from the flrst photometric subsample 
and has a fltted redshift 0.17 ± 0.03 (SDSS-SN internal ID 3077); the second is from the 
second photometric subsample and has a fltted redshift of 0.18 ± 0.04 (SDSS-SN ID 6861). 
Efforts are underway to obtain spectroscopic redshifts for the host galaxies of these events. 
Interpreting the (inflated) photometric redshift errors as Gaussian ( §3.2.11) . the probability 
that at least one of these two candidates has a redshift z < 0.12 is signiflcantly less than 
unity. To be conservative, we assign a systematic uncertainty of +1 SN la based on this 
study. 



•^More pr ecisely, acco r ding t o the photometric typing code, one of the "type A" or "type B" criteria were 



satisfied; see Sako et al 



(|2008l) 
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3.3. Summary of Rate Sample Selection 

In summary, rate-sample selection requirements have been applied to SN candidates 
with z < 0.12 from the 2005 observing season. The resulting sample comprises 16 spectro- 
scopically identified and 1 photometrically identified SNe la. These events are enumerated 
in Table dl and their gri light curves are shown in Figure O along with the best-fit MLCS2k2 
model light-curve. Figure [6] shows the redshift distribution for SNe la for z < 0.21; the lowest- 
redshift photometric candidates with no spectroscopic redshift are in the bin 0.15 < z < 0.18 
and are safely above the redshift cut. 
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Fig. 5. — gri light curves for SNe la used in this rate measurement. Black points show the 
SDSS SN photometry from SMP. The errors on the photometry are shown. Solid curves 
denote the best-fit SN la model light curves in g (green), r (red), and i (dark red) from the 
MLCS2k2 light-curve fitter, and corresponding dotted curves show the 1-sigma model error 
range. The curves and data for the different passbands have been vertically offset for clarity. 
The flux offsets are the same for each SN. 
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z = 0.0570 




-20 -10 10 20 30 40 50 
Observer Frame Epoch [MJD - 53664.5] 

Fig. 5. 
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— Continued. 
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Observer Frame Epoch [MJD - 53686.1] 



Fig. 5. — Continued. 




Fig. 6. — Redshift distribution for events passing the rate- measurement selection require- 
ments in §3.11 Contributions include: spectroscopically confirmed SNe la (black), photo- 
metric SNe la with host galaxy redshifts (gray), and photometric SN la candidates with no 
spectroscopic redshift (light-gray). The arrow shows the redshift cut for this analysis. 
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4. Survey Efficiency 

To convert the number of discovered SNe la into a measurement of tlie SN rate, we must 
liave an estimate of tlie efficiency for discovering SNe la at 2; < 0.12 that satisfy the sample 
selection criteria. We have two tools at our disposal for this estimate: the artificial SN images 
(fakes) that are inserted into the data stream in real time and Monte Carlo simulations of 
the 2005 observing season. 

4.1. Use of Artificial SN Images 

As noted in §2.31 the fake SNe images are used to measure the efficiency of the on- 
mountain software pipeline for point-object detection on a variety of galaxy backgrounds 
and observing conditions. The fake SNe are also used to measure the efficiency of human 
scanners for identifying objects as SNe. While the fakes were designed to model realistic 
SN la light curves, the dependence on the redshift distribution results in only 18 fake 
light curves with redshift z < 0.12. Although all 18 low-redshift fakes were recovered by 
the SN search pipeline, using such a small sample to measure the pipeline efficiency would 
result in large statistical and systematic uncertainties. Furthermore, the fake light curves 
were generated with distributions of Ay and Amis that were not realistic, which complicates 
the interpretation of discovery efficiency as a function of redshift. 

To obtain a more reliable determination of the survey efficiency, we use fake SN la at 
all redshifts in the following way. We first use the fakes to measure the object-detection 
efficiency as a function of the signal-to-noise ratio (SNR) in the (?, r, and i passbands. The 
detection efficiency, defined as the ratio of the number of fake epochs detected as objects 
by the on-mountain software pipeline to the number of fakes inserted into data images 
at a given signal-to-noise, is shown in Fig. 7. While the object detection efficiency as a 
function of magnitude or redshift is sensitive to observing conditions (seeing, clouds, moon), 
the efficiency as a function of SNR is robust against such variations in conditions. As a 
check that the SNR is an adequate parameterization of the point-source detection efficiency, 
we have split the sample of fakes into a low-redshift and a high-redshift subsample and 
determined the efficiency as a function of SNR for each set independently. We find that 
the results are consistent. With the efficiency as a function of signal-to-noise ratio known, 
one can estimate the SN discovery efficiency as a function of redshift for any choices of SN 
la light-curve models, observing conditions, and population distributions. These efficiency 
functions measured with fakes are used in the Monte Carlo simulation ( §4.2p to verify that 
the software pipelines were fully efficient at low redshift. 
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Fig. 7. — The mean object detection efficiency as a function of signal-to-noise ratio for 
SDSS gf-band (top), r-band (middle), and i-band (bottom). The efficiency is derived by 
coimting the fraction of fake images detected by the survey pipeline. The binomial errors on 
the efficiency measurements arc shown. The solid lines show the result of a polynomial fit to 
the efficiency measurements. These efficiency functions are used to simulate the difference 
imaging software in the Monte Carlo simulations of the SN hght curves. 
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The fakes also provide information on the efficiency of the human scanners to correctly 
label as possible SN candidates those fakes that were detected as objects by the software 
pipeline. For the 2005 season, 91% of all epochs of fakes visually scanned by humans were 
flagged as SN candidates, and 95% of all detected fake SNe were flagged by humans as 
SN candidates at least once. The 5% of fakes that were never identifled by humans as SN 
candidates were detected on only a single epoch by the software pipeline, either because 
they were at high redshift or because they reached peak light well before or well after the 
observing season. Essentially all fakes detected on two or more epochs by the software 
pipeline were flagged by humans as SN candidates at least once. Given the selection cuts in 
§3.H the human scanning efficiency is 100% for SNe la contributing to this low-redshift rate 
measurement. 

Summary information on the efficiency of the software pipeline and the human scanners 
to detect fakes is presented in Figure [H which shows the detection efficiencies, i.e., the 
fraction detected by the pipeline and the fraction identifled as SN candidates by humans, 
vs. peak (yf-magnitude. The arrows indicate the peak ^f-band magnitudes for an unextincted 
normal and for an unextincted sub-luminous 1991bg-like SN la at z = 0.12, according to the 
MLCS2k2 model. This flgure indicates that, for the assumed SN la model used to generate 
the distribution of fakes, the combined software+human detection efficiency is essentially 
100% for SNe la in the redshift range z < 0.12. 
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Fig. 8. — The efficiency for identifying fakes as SN candidates, as a function of peak g^-band 
magnitude. The dashed curve is the efficiency for detection by the software pipehne, and 
the dotted curve is the efficiency for evaluation of the fakes as SN candidates by the human 
scanners. The arrows indicate the peak magnitudes for a normal and for a 91bg-like SN la 
at a redshift of 0.12 according to the MLCS2k2 model. The binomial errors on the efficiency 
are shown. 



4.2. Monte Carlo Simulations 



To determine the SN la selection efficiency with high precision and to study systematic 
uncertainties for the rate measurement, we have developed a detailed Monte Carlo light 
curve simulator (MC). The MC simulates individual light curve data points based on real 
observing statistics, but without the added complexity of adding fake SNe to images. The 
MC light curves can be generated and analyzed much more rapidly than the fakes, so the 
MC can be used to rapidly simulate very large numbers of SN la light curves to estimate the 
SN discovery efficiency and the uncertainty in the efficiency due to assumptions about the 
SN la model distributions. The MC code uses the MLCS2k2 model to generate simulated 
SN la light curves instead of the stretch/ Amis model that was used to generate the fakes. 

For each simulated SN la, the following parameters are randomly drawn from parent 
distributions: 



redshift, z: 

Drawn from a distribution proportional to the comoving volume element, which as- 
sumes a constant SN la rate per unit comoving volume. 

host galaxy extinction, Ay: 

Drawn from a distribution P{Ay) oc e^"^^/^, with r = 0.4. The ICardelli et al.l ( 19891 ) 



reddening law, with Ry = 3.1, is used to extrapolate the e xtinction t o othe r wave 



lengths. The choice of r = 0.4 was guided by the studies of iJha et al.l (120071 ) and is 
consistent with the inferred extinction distribution for spectroscopically conffimed SNe 
la in the SDSS SN sample. As we discuss later in this section, the exact choice of r 
makes no practical difference to this rate measurement. 

3. MLCS2k2 light curve shape/luminosity parameter, A: 

Drawn from a bimodal Gaussian with a standard deviation of 0.26 for A < and 
0.12 for A > 0, and truncated to lie within the valid range of the MLCS2k2 model, 
—0.35 < A < 1.8. The bimodal Gaussian is based on study of the conffimed SNe la 
in the SDSS-II Supernova Survey ffist year data. 

4. time of peak light in rest-frame 5-band: 

Drawn randomly from the interval 53616 < MJD < 53705 (Sept. 1 - Dec. 1, 2005). 

5. sky position: 

Drawn randomly from the range of the survey. 



6. 



location within host galaxy: 

Drawn from a distribution proportional to the host galaxy surface brightness (see 
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below). This variable is used only to determine galaxy background light, not extinction. 

We note that simulated photometry is generated only at epochs for which we obtained 
photometric imaging at the corresponding sky location, and therefore the determination of 
the selection efficiency naturally accounts for the temporal inhomogeneity in sky coverage, 
as can be seen in Figure [TJ 

Using these parameters for each SN, rest-frame UBVRI magnitudes are generated from 
the MLCS2k2 model for all dates on which the survey took data at the selected sky position. 
These magnitudes are modified according to the host galaxy extinction, i^-corrected to the 
observed SDSS gri passban ds, and further mod ified according to the estimated Milky Way 



extinction at that position (jSchlegel et al.lll998l ). The zero-points from the survey are used 



to convert the gri magnitudes into flux values that would have been measured in ADUs 
by the SDSS 2.5 m telescope. The CCD gains are then used to determine the number 
of photo-electrons, and hence the signal and noise. Additional noise is computed for each 
measurement based on the measured observing conditions at each epoch, in each passband, 
at the assigned sky location. Sky noise is simulated by integrating the estimated sky noise 
per pixel over an effective aperture with a size determined by the local PSF estimate from 
PHOTO. Noise from the host galaxy is sim ulated by associatin g the SN with a host from the 



SDSS galaxy photometric redshift catalog (lOyaizu et al.ll2007l ) sele cted to have a photometric 



redshift equal to the assigned SN redshift. In the SD SS DR5 (lAdelman-McCarthy et al. 



20071 ) photoPrimary database (jStoughton et al.ll2002l ). each such galaxy image is flt with 



both an exponential and a de Vaucouleurs surface brightness profile. We use the exponential 
model in the r-band as a probability distribution from which the SN position within the 
galaxy is drawn. That is, the galaxy noise model assumes that the SN la rate is roughly 
proportional to r-band stellar luminosity. The estimated contribution of the galaxy light to 
the noise in each passband is computed by convolving the exponential galaxy model with 
the PSF in the survey image. In practice, this procedure is computationally expensive, so we 
pre-compute the noise values on a grid of model parameters and perform a multi-dimensional 
linear interpolation to obtain an estimate of the galaxy noise. 

As a consistency check of the MC as a representation of the SN data, we compare 
the distributions of signal and noise in gri for the MC sample to the signal and noise for 
all photometric epochs for the low-redshift SNe la in the rate-measurement sample. The 
comparison of the distributions is shown in Fig. O The distributions of signal and noise are 
in good agreement, indicating that the MC model, and the assumed parameter distributions 
therein, provide a reasonable representation of the low-redshift SN la sample. 
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Fig. 9. — For SN la that satisfy the selection requirements the signal-flux and noise 
distributions are shown for all photometric epochs in the g,r,i filters. Each distribution is 
shown for SDSS data (dots) and for the simulation (histogram) that has been scaled to have 
the same number of entries as the data. 
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Having shown that the MC generates photometry that is consistent with our observed 
low-redshift SN la sample, we can use the MC to provide a more rehable determination of 
the detection efficiency of the on-mountain search pipehne than we could obtain with the 
relatively small number of low-redshift fakes. For each epoch generated by the MC, we use 
the efficiency as a function of signal-to-noise ratio in each passband (from Fig. [7j) to determine 
the efficiency of the search pipeline to detect SNe la at all redshifts. The resulting software 
pipeline detection efficiency as a function of redshift, based on a MC study using 15640 
total SNe (920 in each of 17 redshift bins) is shown in Figure [TUl the efficiency is 100% over 
the redshift range z < 0.12. Since the fakes tests of §4.11 showed that the human scanning 
process causes a negligible loss of efficiency, we conclude that the combined efficiency for 
SN detection by the pipehne and identification as a candidate by humans is 100% over the 
redshift range of interest. This does not guarantee that the efficiency of the photometric 
typing code used for spectroscopic target selection ( §2.21) is also 100%, but the studies of 
§3.21 indicate that, with the exception of SN 9266, there were no losses due to the target 
selection algorithm. 

The final step is to use the MC to compute the survey discovery efficiency e for a SN la 
sample defined by the selection requirements in §3.11 This efficiency is the ratio of the 
number of SNe la that are detected by the pipeline, identified by humans, and that pass 
the selection criteria of §3.11 to the total number of SNe la that reach peak light during the 
survey, i.e., between Sept. 1 and Nov. 30. While we have seen that the detection efficiency 
is essentially 100% out to 2; = 0.12, e is less than 100% primarily because the selection 
requirements on light curve coverage (cuts 4 and 5 in §3.ip remove some SNe la that peaked 
in early September or late November. Using the MC sample of 15640 light-curves mentioned 
above and fitting a linear function to the resulting selection efficiency in the redshift range 
< 2; < 0.12 gives 

e{z) = (0.78 ± 0.01) + (-0.13 ± 0.14) z (1) 

That is, the survey efficiency is approximately constant at low redshifts, changing by only 
~ 1% over the redshift range of the rate measurement. The mean SN la discovery efficiency 
for our rate sample is (e) = 0.77 ± 0.01. 

While the data-MC comparison in Fig. [Hlindicates that we have made a consistent choice 
of the parameter distributions for the MC model, to estimate the systematic uncertainty in 
the discovery efficiency we vary the assumed MC parameter distributions and recompute the 
efficiency. We find that varying r, the parameter controlling the extinction distribution, has 
the largest systematic effect on the determination of the discovery efficiency from the MC. 
Varying r over the range 0.2—0.6, the estimated discovery efficiency for the rate-measurement 
SN la sample changes by less than a percent. 
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Fig. 10. — Search pipeline software efficiency as a function of SN la redshift, as determined 
from the Monte Carlo simulation. 
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5. SDSS SN Results 
5.1. Volumetric SN la rate 



For the purpose of interpreting the SN observations as a volumetric SN rate, we will 
assume a spatially fiat cosmology with non-relativistic matter density Qm = 0.3, dark energy 
density = 0.7, and dark energy equation of state parameter w = p/ p = —1. For the 
low-redshift rate measurement presented in this section, the dependence on cosmological 
parameters of the survey selection efficiency is negligible, and so the uncertainty in the rate 
due to uncertainty in cosmology is due entirely to the difference in the volume of the survey. 
A change in Qm of 0.02 would lead to a 4% change in the rate measurement. 

The observed volumetric SN rate, ry, is defined as 



ry = ^ , 2 
VTe 

where N is the number of SNe in the sample, and VTe is the effective product of the survey 
volume, V, the observer- frame survey duration, T, and the SN discovery efficiency, e{z), 
estimated in §li2|, 

rZmax J^. 1 

VTe = (0T) / dz 7^ • (3) 

Here B is the solid angle covered by the survey and u{z) is the comoving distance in the 
Friedmann- Robertson- Walker metric, 

io H{z') Ho Jo ^nUl + z'^ + ^A 

If the survey efficiency is independent of redshift, and if the redshift range covered by the 
SN observations is small, then VTe ~ (VTe) /{I + (z)), where V is the survey volume and 
{z) is the volume-weighted mean redshift of the survey. 

For the SDSS-II Supernova Survey we have N = 17, Zmin = 0, z^ax = 0.12, {e{z)) = 
0.77 ± 0.01, T = 89 days = 0.244 years, and 9 = 0.08277 * 0.98 steradians. This value for 6 
is 98% of the actual sky area covered by the survey, due to the masking of bright stars and 
variable sources. Substituting these values into Eqn. ([2]), we find a volumetric SN la rate of 
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r-v = [2.93l°:J^(systematic)l°:?'i' (statistical)] x 10~^ SNe Mpc"^ /i^g year"^ , (5) 

with hjQ = Hq (70 km s~^Mpc~^)^^ and Hq the present value of the Hubble parameter. The 
statistical errors quoted represent the standard frequentist 68.27% central confidence interval 
on the mean of a Poisson distribution. The systematic uncertainty represents uncertainty on 
our determination of the SN selection efficiency ( §4.21) and on the number of photometrically 
identified SNe la ( §3.21) . This measurement represents the volume-averaged SN la rate at 
z < 0.12. When rate measurements are plotted vs. redshift, it is generally assumed that 
the rate is constant over the sampled redshift interval. If we assume that the SN la rate is 
constant at z < 0.12, then Eqn.([5]) can be interpreted as the rate at the volume-weighted 
mean of our redshift range, (z) = 0.09, and we make this assumption when plotting the 
result. Our result is shown along with previously reported SN la rate measurements in 
Fig- [m but we defer discussion of comparison and combination with other measurements to 



5.2. SN la Rate per unit galaxy luminosity 



Early measurements of the SN rate were generally derived from SN observations that 
targeted known galaxies; for these surveys, the SN rate is most nat urally measu r ed as a rate 
per unit luminosity in some passband, traditionally the 5-band. iBlanc et al.l (120041 ) have 
converted a number of measurements from the literature of the SN l a rate per unit lun i inosit y 
to rates per unit yolurn e, and in Tableljlwe ad opt their values for the lCappellaro et al.l (119991 ): 



Madewick et al.l (120031 ): iHardin et al. 



(120001 ) rate measurements. 



For completeness, we convert our volumetric rate to a rate per unit galaxy luminosity in 
th e SPSS passba n ds. T he galaxy luminosity functions in the SDSS passbands are estimated 



m 



Blanton et al.l (120031 ). The corresponding luminosity densities in the ugriz passbands, 
at a mean redshift of {z) = 0.1 are 1.60 ± 0.32, 1.25 ± 0.05, 1.29 ± 0.04, 1.48 ± 0.05, and 
1.89±0.05, in units of lO^LQ/iyoMpc"^, where Lq is the solar luminosity. In combination with 
the volumetric rate measurement of Eqn.dS]), this yields the SN la rate per unit luminosity 
in the SDSS passbands, {rL)ugriz/h^7o = O-'^^Stom SNu„, 0.235tofe SNu^, 0.227i[}:|]^ SNu^, 
O.mtof^ SNui, and O.lbQtoM SNu„ where 1 SNu^ = 1 SN 10"^°^^ (100 yr) 
the luminosity in passband x, in units of solar luminosities. 



^, with Lg 
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5.3. SN la Rate as a function of host galaxy type 



Recent measurements have shown that the specific SN la ra t e is hi gher in star-forming 
galaxies than in passive galaxies. For example, iMannucci et al.l (120051 ) found that the SN 
la rate p er unit stellar mass is ~ 20 — 30 times higher in late-type galaxies than in E/SO 
galaxies. ISuUivan et al.l (120061 ) have found a similar trend in the SNLS data. We will consider 
the trend of SN rate vs. star formation activity using the SDSS-II Supernova Survey sample 
in a forthcoming publication. 

Here, we consider the low-redshift SN la rate vs. host galaxy type. We have consid- 
ered sever al photometric galaxy-type indi cators that are accessibl e through the SPSS DR5 
database (lAdelman-McCarthy et al.l 120071 ). including u — r color (IStrateva et al.ll200ll ): the 
likelihood of the de Vaucouleurs model fit to the galaxy surface brightne ss profile relative 



to that of an exponential n iodel fit; and the (inverse) concentration index (jShimasaku et al. 



2001: Yamauchi et al. 



2005), defined as the ratio of the radii that contain 50% and 90% of the 



Petrosian fiux (see IStoughton et al.l (120021 ) for definitions of these quantities). These param- 



eters are listed in Table [3] for the host galaxies of the SNe included in the rate-measurement 
sample. Note that the u — r color in Table [3] is computed using SDSS model magnitudes 
( IStoughton et al.l |2002|) and is not i^-corrected to the galaxy restframe. A host galaxy is 
associated with each SN by determining the nearest object, based on a measure of the galaxy 
image size. Specifically, the SDSS DR5 catalog includes the parameters of an iso-photal el- 
lipse for each galaxy-like object, and for each of the ugriz filter bands, defined as the e llipse 
where the object surface brightness is 25 magnitudes arcsec"^ (IStoughton et al.ll2002l ). We 
define the distance to a potential host galaxy to be the semi-major axis of the ellipse that 
is similar (has the same aspect ratio and orientation) to the r-band iso-photal ellipse and 
that intersects with the position of the SN. The host galaxy for each SN is defined to be the 
nearest object in this measure. For the SNe listed in Table [H the association of SN with 
host galaxy was confirmed through visual inspection of the images. 

Early-type galaxies generally display red colors (large u — r), are reasonably well fit 
by a de Vaucouleurs surface brightness profile (large values of the relative de Vaucouleurs 
likelihood), and show relatively strong central light concentration (low values of the inverse 
concentration index). Consequently, these three indicators are strongly correlated and tend 
to give a consistent classification into early and late photometric types, indicated by the last 
column in Table [31 However, the classifications based on the three indicators do not always 
agree, in which case we have made a judgement based on visual inspection of the galaxy 
image and, where available, a high signal-to-noise galaxy spectrum. For these five cases, the 
host type indicated in the last column is marked with an asterisk. 
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Table 3. Host Galaxies for SNe la in the rate sample. 



SDSS Id 


(T ( T2000 0"! 


A ( T2000 0"! 


Uj I 


n P 1 1 POl 1 1 Pll Tc; 
1 1 kpIi noon 


POD PPTli'T';^ "hlOTl 
111 (iPV 


rTocif" T^A/'"nP 
J-J-VJOu J- VUC 


1241 


22 30 41.15 


-00 46 34.5 


2.82 


0.91 


0.382 


Early 


1371 


23 17 29.70 


+00 25 46.8 


2.97 


0.00 


0.377 


Early* 


2561 


03 05 22.64 


+00 51 35.0 


2.59 


0.00 


0.410 


Late* 


3256 


21 57 04.19 


-00 13 24.5 


1.99 


0.00 


0.442 


Late 


3592 


01 16 12.71 


+00 47 26.0 


2.21 


0.00 


0.427 


Late 


3901 


00 59 24.11 


+00 00 09.5 


1.40 


1.00 


0.418 


Late* 


5395 


03 18 33.80 


+00 07 24.0 


1.29 


0.49 


0.360 


Late* 


5549 


00 12 59.97 


+00 14 54.9 


1.01 


0.46 


0.461 


Late 


5944 


01 56 48.50 


-00 12 45.3 


2.57 


0.00 


0.469 


Early* 


6057 


03 30 12.89 


-00 58 28.1 


1.79 


0.00 


0.485 


Late 


6295 


01 34 41.84 


-00 36 15.2 


2.97 


1.00 


0.312 


Early 


6558 


01 26 48.46 


-01 14 17.3 


2.23 


0.00 


0.427 


Late 


6962 


02 35 26.58 


+01 04 28.3 


2.71 


1.00 


0.379 


Early 


7147 


23 20 04.44 


-00 03 20.2 


3.22 


1.00 


0.350 


Early 


7876 


01 16 43.87 


+00 47 36.9 


1.71 


0.00 


0.532 


Late 


8719 


00 30 53.23 


-00 43 07.3 


1.12 


0.01 


0.412 


Late 


9266 


03 20 43.19 


-01 00 08.2 


2.25 


0.00 


0.398 


Late 



*At least one of the three photometric type indicators indicates a different type from that 
listed. 



Note. — SDSS Id denotes internal candidate designation, a and 5 are the coordinates of 
the host galaxy of the SN. The photometric morphology indicators, u — r, de Vaucouleurs 
likelihood, and concentration index are described in ( §5.31) . 
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Of the three photometric type indicators, u — r correlates most strongly with the host 
type we have assigned to each galaxy in the las t column of Table p The distribution of 
SDSS galaxies is approximately bimodal in u — r (IStrateva et al.ll200ll ). suggesting a natural 
division between early and late types. We therefore use m — r as the galaxy classifier for the 
purpose of studying the SN rate vs. galaxy type. This is preferable to using the 'host type' 
classification in Table [3l since the subjective human judgem ent required to determine the 
latter makes it difficult to determine its population properties. IStrateva et al.l (1200 ll ) suggest 
that u — r = 2.2 is an optimal separator between e arly and late types. However, our catalog 
of galaxies with photometric redshifts in stripe 82 (lOyaizu et al.l 120071 ) appears to be better 
separated into two subpopulations using u — r = 2.4. Since a division at m — r = 2.4 also 
provides better agreement with the subjective 'host type' classification in Table [3l we use 
this color cut to separate the hosts into early [u — r > 2.4) and late {u — r < 2.4) types for the 
relative rate measurement. Using a large sample of galaxies from the SDSS DR5 database, 
we find that the fractional r-band luminosity densities for early and late-type galaxies at 
redshifts ^ < 0.12 are 54% and 46%, respectively. From Table [3l we find that the SN la rate 
per unit r-band luminosity is ~ 1.68]'^o 4^ tir nes higher in la t e- type galaxies than in early-type 
galaxies. Using the luminosity functions of iBlanton et al.l (120031 ). we find that the absolute 
rates per unit luminosity are ri/hjQ = 0.085]'^o;q2 SNu,. (early) and 0.142lQ;g3 SNu^ (late). 
The evidence for a larger SN la rate in late-type galaxies is statistically marginal with the 
current low-redshift sample. The systematic uncertainty is also significant: if we place the 
host galaxy type cut at m — r = 2.2, we find no significant difference between the rate per 
unit luminosity in early- and late-type hosts. 



6. Fitting SN Rate evolution models 

As noted in ^ models for SN la progenitors in principle can be distinguished by their 
predictions for the evolution of the SN la rate with cosmic time. In this section, we present 
a general maximum likelihood method of fitting SN observations to models with a redshift- 
dependent SN rate. We then apply the method to a recently discussed SN la rate model, 
using data from the SDSS-II Supernova Survey and from other published rate measurements. 



6.1. Maximum Likelihood method 



In this section we describe a method for fitting SN data to models of the SN rate withou t 
bin ning the data. Th e meth od is similar to the methods described in IStrolger et al.l (120041 ). 
and IStrolger fc Riesd (120061 ). A distinguishing feature of our analysis is that it allows for 
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combining multiple data sets and accounts for systematic errors. 

A general model for the volumetric SN rate can be written as rv{z; p), where p represents 
the set of model parameters. According to the model, the total number of detected SNe 
follows a Poisson distribution with mean value 



(7V(p))^^~d.ere(.)!l(nP),.^w|, (e) 

where all symbols were defined in ^ The probability of detecting a SN at redshift z is given 
by the integrand of Eqn.Q, P{zi) oc d{N{p)) /dz, giving a likelihood function for detecting 
SNe at the N observed redshifts {zi}, 



_ e-(^(p))(Ar(p))^ - I d{Nip)) 

^^^^^^'P^- m l\{N{p)) dz ■ 

The corresponding log-likelihood function, suppressing terms that do not depend on the 
model parameters, is 



logL({;.,};p) = -(iV(p)) + f^log ( eTe(z.) '-^^^ 



U [Zi 



du 
dz 



The best-fit model is determined by maximizing the log-likelihood with respect to the model 
parameters, p. To incorporate information about the systematic error in our fits, we weight 
the contribution to the log-likelihood for each data set by multiplying each term in the log- 
likelihood function by the factor cr'^tat/i^stat + ^lyst)^ where astat and asyst are the statistical 
and systematic errors for the measurement. This factor assumes that the systematic errors 
are approximately Gaussian and independent of the statistical errors. We note that auxiliary 
information about the model parameters and uncertainties in the survey parameters B, T, 
and e{z) could be incorporated in a more rigorous way via prior probability distributions. 
However, this would require full knowledge of the probability distribution functions for the 
efficiency, subject to all possible variations of systematic effects, which is in practice unknown. 

To combine data from multiple surveys, the log-likelihood functions for each survey are 
added together, using the appropriate values of e{z), 9, and T for each survey. The advantage 
of this method is that it does not involve binning the SN data in redshift; however, it does 
require knowledge of the efficiency function e{z) for each survey. To evaluate the goodness of 
fit of a given model, one can use, e.g., the Kolmogorov-Smirnov (KS) test applied to the data 
and to a large-statistics Monte Carlo sample generated from the best-fit model parameters. 
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The code for fitting rate models was tested on large MC samples (~ 1000 SNe), and the MC 
model parameters were accurately recovered. 

As an illustration of the likelihood method, we apply it to the SDSS-II Supernova 
Survey data, assuming a redshift-independent model over the redshift range probed by the 
data, rv{z) = constant. In this case, the rate that maximizes the likelihood can be shown 
analytically to be given by Eqn. [2l The probability for this model from the KS test statistic 
is Pks = 0.42, meaning that if the model is correct, 42% of sample observations drawn from 
the model would have a KS test statistic as large or larger than that found in comparing this 
data set to the model. In the discussion that follows the probabilities from the KS test are 
given as rough estimators of the goodness of fit only; the distribution of the KS test statistic 
does not in general have an analytic form when model parameters are estimated from the 
data. 



6.2. SN Rate Models and Star Formation History 

As discussed in ^ measurements of the SN la rate provide a means to distinguish 
between models of SN la progenitor systems. The connection between the observed SN la 
rate and the progenitor systems is made through the relation of the SN rate to the cosmic 
star formation history. Sometime after a population of stars form, a fraction of them will 
end up in binary systems that are producing SN la explosions. If we denote the distribution 
of delay times between formation of the progenitor systems and the SN explosions by D(t), 
then the volumetric SN la rate rv{t) and the cosmic star formation rate p{t) are related by 



Mt) 



dt'p{t')D{t - 1') 



(9) 



We can therefore constrain models for the distribution of delay times, D{t), by comparing the 
SN la rate and the star formation rate. A discu ssion of predict ed delay-time distributions for 
a variety of SN la progenitor models is given in lGreggid ( 20051 ). A simple model distribution 
that allows for two distinct contributions to the SN la rate is 



D{t) =A + B6{t) 



(10) 



where 5{t) i s the Dirac delta functio n . Thi s 'A + 5' model was proposed by lMannucci et al. 
(120061 ) and IScannapieco &: BildstenI (I20051) and it has been used in SN rate studies by the 



SNLS (INeill et al.l 120061 ) and 



Sullivan et al.l (120061 ). The SN rate can be written ry(t) 
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Ap(t) + Bp{t), where p{t) is the stellar mass density. The B term represents an instantaneous 
or prompt SN la component and the A term represents an extended component in which SNe 
la form with uniform probability in the time interval following star formation. In addition 
to the 'A + 5' model, we also consider a simple model in which rv{t) evolves as a power law 
in redshift, independent of considerations of star formation history. 



6.3. Rate Measurements: Combining Data Sets 

The constraints on redshift-dependent models of the SN la rate are improved if one uses 
SN observations over a wide range of redshifts. In the following, we combine the low-redshift 
rate measurement from the SDSS-II Supernova Survey with other SN la rate measurements 
in the literature. For each data set, we require both the SN redshifts and an estimate of the 
redshift-dependent selection function e{z), and we therefore restrict ourselves to using data 
sets for which it is straightforward to infer the redshift depe ndence of the eff i ciency . We 



note that several authors, including iBarris fc Tonryl (120061 ) and iPoznanski et al.l (120071 ) have 



made SN rate measurements based on samples of photometrically identified SNe. However, in 
combining data sets for the present analysis, we will restrict ourselves to rate measurements 
that are based primarily on spectroscopically identified SNe. Of the nine previously published 
rate measurements that have been based on primarily spectroscopically identified SNe, shown 
in Fig. [m we will make use of four, in addition to the one in this work. These five rate 
measurements are shown in bold font in Table |H The weighting factors, used to account for 
the systematic uncertainty on each measurement, are listed in the last column of Table |H 
In cases where the uncertainty on the measurement is asymmetric, we define the weighting 
factor to be the mean of the upper and lower weighting factors. Varying the weighting factor 
between the extremes of using the smaller weight, and of using the larger weight, the best-fit 
parameters change by ~ 5% of the statistical error. In the subsections below, we briefly 
describe the data from other measurements that we include in the model fits and how we 
describe their efficiency function. We also discuss measurements that we exclude from the 
model fits. 
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Table 4. SN la Rate Measurements. 



Reference 



Redshift 



Mean 



SNe 



Rate 



stall '-'tot 





Range 


Redshift 




[10-5 SNe /ifo Mpc-3 yr-i] 




CaoDellaro et al. (^1999)* 


~ 


~ 


70 


2.8 ±0.9 


N/A 


This work* 


0-0.12 


0.09 


17 




0.988 


Madffwick et al. ('2003^ 


0-0.19 


0.10 


19 


3.1 ± 1.6 


N/A 


Blanc et al. ('2004) 


0-0.3 


0.13 


14 


9 n+0-84 
^•^-0.72 


N/A 


Hardin et al. C2000) 


~ 0.02 - 0.2 


0.14 


4 


3.4t?:? 


N/A 


Daiilen et al. C2004) 


0.2 - 0.6 


0.45 


3 




N/A 


Neill et al. (2006)* 


0.2 - 0.6 


0.45 


73 


4 2+^-'' 
^•^-1.1 


0.492 


Tonrv et al. (2003) 


~ 0.25 - 0.6 


0.46 


8 


4.8 ± 1.7 


N/A 


Pain et al. (2002^" 


0.25 - 0.85 


0.55 


37 




0.643 


Dahlen et al. (2004^ 


1.0 - 1.4 


1.2 


6 


n.5ttl 


0.686 



The value of the rate has been corrected to our assumed cosmology, according to equation 3 of lPain et al 



(12002). 



Note. — Measurements included in the model fits are shown in bold face, and are marked with an 
asterisk. Mean redshift refer s to the mean of t he exp ected SN redshift distribution, under the assumption 
of a constant SN rate. For iMadgwick et al.l (|2003l l. this is estimated as the mean of the observed SN 
redshift distribution. Systematic and statistical errors, when reported separately, have been combined in 
quadrature. Rate measurements reported here assume constant volumetric rate over the range of each 
survey, agtat is the reported statistical error on the measurement, atot is the sum in quadrature of the 
reported statistical and systematic errors. 
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Fig. 11. — Measurements of the SN la rate discussed in §6.31 The SDSS-II Supernova 
Survey measurement in this paper is shown as the sohd black square. Measurements for which 
the data is used in the model fits are shown as solid circles (see Table Hj), and measurements 
not used in the fits as open circles. To plot each measurement, we have assumed in each case 
a model in which the rate is constant over the redshift range covered by that measurement. 
The rate as a function of redshift for the best fitting M + 5' and power law models are 
overlaid. 
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0.2 0.4 0.6 0.8 



1.2 1.4 
Redshift 



Fig. 12. — Comparison of the observed distribution of SNe and the predicted distributions for 
the A + B and power law rate models. In each panel the shaded region shows the predicted 
redshift distribution of the best-fit model. The figures include the five highhghted data sets 
in Table H 
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6.3.1. Other rate measurement data included in the model fits 



Neill et al.l (120061 ) measured the SN la rate using 73 SNe la from the SNLS. They 
state that their sample is spectroscopically complete, i.e., that e{z) is constant, to z = 0.6. 
Including the solid-angle and survey observation time, the factor QTe{z) = 7.37 x 10~^ 
steradian year. 



The measurement of iPain et al.l (120021 ) is based on data from t h e Sup ernova Cosmology 
Project, covering about 12 square degrees. Although iPain et al.l (120021 ) do not give e{z) 
explicitly, they do provide the redshift distribution of SNe recovered from their Monte Carlo 
simulations which assumed a constant rate per unit comoving volume. With this information, 
we can compute the relative number of MC-generated SNe in each redshift bin and thereby 
the redshift dependence of their efficiency function. Fitting a quadratic function to this 
tabulated efficiency function in the range 0.25 < z < 0.85 gives QTe{z) = (2.68 + 0.61z — 
4:.22z^) X 10~^ steradian year. 



Cappellaro et al.l (119991 ) measured the SN la rate for the local Universe by combining 



data from a number of surveys, including visual searches of nearby galaxies. Although they 
do not provide an efficiency function or a redshift distribution, the redshift range covered by 
the measur ement is so small that we take the quoted result to be the SN la rate at ^ = 0. We 
include the ICappellaro et al.l (119991 ) rate measurement by adding a standard term to the 



log-likelihood function, i.e., a term of the form (ryfO ; p) — ry^cappY I'^a'nr.r^, where ry(0;p) 



is the model prediction at redshift zero, ry^capp is the ICappellaro et al.l (119991 ) measurement 
in Table m and acapp is the quoted error on the measurement. 



Dahlen et al.l (120041 ) measured the SN la rate to z ~ 1.6 using data from the Great 



Observatories Origins Deep Survey (GOODS) carried out with the Advanced Camera for 
Surveys on the Hubble Space Telescope (HST). We use their measurement in the redshift 
range 1.0 < z < 1.4 . Using their scaled efficiency function, as inferred from Figure 14 of 
Strolger et al.l (120041 ). we fit a function A + Bz + Cz"^ + Dz^, valid in the redshift range 
1.0 < z < 2.0. The best-fit parameters are A = — 7.557, S = 55.93, C = —51.07 and 
D = 12.5. This function is normalized so that the number of expected SNe for their survey 
is equal to six, the number they observed in the redshift interval 1.0 < z < 1.4, giving a 
value of &Te{z) = [-10.35 + 76.61^ - 69.95^^ + 17.122^] ^ iQ-^ steradian year. 

The redshift dependence of the efficiency function for the present data set is discussed in 
§4.2( including the solid-angle and survey observation time, QTe{z) = [1.54 — 0.0252;] x 10~^ 



''Just before we submitted this paper, iKuznetsova et al.l (|2007l ) releas ed new SN la rate m easurements 
based on analysis of 57 SNe from HST, including the 42 SNe analyzed bv iDahlen et al.l (|2004l ). 
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steradian year. 



6.3.2. Rate measurement data not included in the fits 



In fitting the models, we ch oose not to include several of the SN la rate measurements 
listed in Table |H We exclude the lDahlen et al.l (120041 ) r ate measurement in t he redshift range 
0.2 < z < 0.6 because the efficiency function, given bv IStrolger et al.l (l2004f). is onlv plotted 
for redshifts greater than 1, and because the 73 SNe from SNLS dNeill et al.ll2006[) in the same 
redshift range dominate the fit in compa rison to the 3 SNe f rom lDahlen et al.l (120041 ). Similar 
reasoning holds for the measurement of iTonry et al.l (120031 ). which is based on 8 SNe in the 
redshift r ange well covered by the SNL S, and for which th e redshifts are not explicitly stated. 
Both the iHardin et al.l (120001 ) and the iBlanc et al.l (120041 ) rate measurements, based on data 
from the EROS microlensing survey, included a requirement that each SN be associated with 
a host galaxy with apparent magnitude R < 19, which introduces a bias against faint hosts. If 
SNe la occurred at a constant rate per unit i?-band luminosity in all galaxies, this would not 
be an issue. However, as noted above, it has been shown that the SN la rat e per unit stellar 



mass (for which the total -R-band lumin osity is a proxy) i s a fun ction of SFR (iMannucci et al. 



20051 : ISuUivan et al.ll2006l ). Finally, the lMadgwick et al.l (120031 ) measurement is based on SNe 
discovered via principal component analysis in spectra obtained by the SDSS galaxy redshift 
survey. This measurement has significant systematic uncertainties that are different from 
those in photometric surveys. In particular, the SNe discovered by this technique must 
lie within approximately 1.5" of the cores of their host galaxies, the radius of the SDSS 
spectroscopic fibers. To derive a SN rate from these spectroscopic observations, assumptions 
are needed about how SNe are distributed within their host galaxies at larger galactocentric 
distances. 



6.4. Fits to SN la Rate Models 



We now consider fits of the combined SN la rate measurements to the rate models 
discussed in §6.21 using the maximum likelihood approach of §6.11 The errors quoted below 
are the values of the fit parameters for which the log-likelihood function changes by 1/2 
compared to its maximum, which assumes that th e likelihood function is approximately 
Gaussian. We use the MINUIT software package (jjames &: Rood |l99J) for the function 
optimization and error analysis. 
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6.4.1. Power-law redshift evolution of ry 

We first consider a simple two-parameter model that describes power-law redshift evo- 
lution of the SN rate independent of consideration of the star formation history, rv{z) = 
a{l + z)^. The best-fit power-law model is shown as the dashed curve in Fig. [11], and the 
predicted redshift distribution is shown in Fig. [T21 Fitting this model to the five data sets, 
we find 

a = {2.Qtoi) X 10"^ SNe Mpc"^ h^o yr~^ 
p = (1.5 ±0.6) 
= -0.80 

where p^/? is the correlation coefficient between the two fitted parameters. The KS probability 
for this model is pxs = 0.63. We emphasize that the fitted value of /5 is greater than 0, i.e., 
the rate is determined to be an increasing function of redshift, at the ~ 2.5a level. 



6.4.2. The 'A+B' model 



We next consider the 'A + B' model, with D(t) given by Eqn. [TUJ As discussed 
by lForster et al.l (120061 ) . there is still significant uncertainty on the cosmic star formation rate 
(SFR), which is a limitation for placing observational constraints on SN delay time models. 
In what follows we choose one estimate of the SFR, and do not propag ate the systematic 
uncertainties in the SFR. We follow the appr oach of iNeill et al.l (120061 ) and take the star 
formation rate from lHopkins Sz BeacomI (120061 ). The functional form of the star formation 
rate is 



a + bz 
1 + {z/cY 



hoo Mq yr-i Mpc- 



where hwo = Hq (100 km sec'^ Mpc"')-\ a = 0.0118,6 = 0.08, c = 3.3, and d = 5.2. For 
the stellar mass density (the A com ponent) we integrate the star formation rate over time; 
as mentioned by iNeill et al. J2OO6I ). this can be expected to overestimate the total stellar 
mass density relative to estimates of the stellar mass density that are based on luminosity, 
as it includes a contribution from stars that have burned out. Performing the fit using the 
five data sets gives 



A = (2.8 ± 1.2) X 10 

B = (9.3_3.i 
PAB = -0.78 



SNe M"^ yr"^ 



X 10"^ SNe M"^ 
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where pab is the correlation coefficient between the two fitted parameters. The KS prob- 
abihty for this model is pks = 0.71. The best-fit 'A + 5' model is shown as the solid 
curve in Fig. [TT], and the predicted redshift distribution is shown in Fig. [T21 We note 
that the uncertainties on the A and B parameters are ~ 43% and ~ 35%, respectively. 
For comparison, if we perform a fit to the 'A + 5' model, suppressing the SDSS data, the 
uncertainties on the fit parame ters are ^ 53% an d ~ 38%, respectively. Our analysis here 
is similar to that presented by iNeill et al.l (120061 ). with the primary differences being that 
we use a different subset of the available data, and we u se a maxirnum lik elihood method 
to fit the data to models of the SN rate. For comparison, iNeill et al.l (120061 ) found values of 
A = (1.4 ± 1.0) X 10"^^ SNe M^^ yr'^ and B = (8.0 ± 2.6) x 10"^ SNe Mq\ Both analyses 
find evidence for two components to the SN rate with the significance of the 'A' (extended) 
component less than that of the '5' (prompt) component. 

We note that one cannot accurately judge the goodness of fit of this model using a 
visual inspection or fit to Fig. [TTl since the measurements are each plotted assuming 
a constant-rate model. A better picture of the goodness of fit is given by Fig. [121 which 
shows the observed redshift distribution for the five data sets compared with the predicted 
redshift distributions for the '^4-1-5' and power-law rate models convolved with the measured 
efficiency functions for the different measurements. The agreement between the predicted 
distributions for both evolving models and that of the data is quite reasonable. 



7. Conclusions 

We have presented a measurement of the SN la rate in the redshift range < z < 0.12 
from the first season of the SDSS-II Supernova Survey. After selection cuts, the rate- 
measurement sample includes a total of 17 SNe la, of which 16 were spectroscopically con- 
firmed. The final SN in the sample is a highly extincted, photometrically identified SN la 
with a measured host galaxy redshift. The insertion of artificial SNe in the data stream and 
the use of detailed Monte Carlo simulations of the survey efficiency, along with the rolling 
nature of the SDSS-II Supernova Survey, have enabled us to obtain a SN la rate measurement 
with smaller systematic uncertainties than previous measurements in a comparable redshift 
range. 

We have also applied a maximum-likelihood technique, which enables us to account for 
systematic errors and to fit multiple SN data sets to models of the SN rate as a function 
of redshift. This maximum likelihood method makes optimal use of the available data, but 
requires estimates of the SN detection efficiency, and its uncertainty, as a function of redshift. 
We have applied this technique to a combination of recent SN la data sets, focusing on the 
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'y4 + 5' model that relates the SN la rate to the cosmic star formation rate. 

Models in which the SN la rate evolves with redshift are preferred over a model with 
a constant rate, but the data do not distinguish significantly between a simple power-law 
evolution of the SN la rate with redshift and the 'A + B' model. The A an d B parameter 



values we obtain are in good agreement with the results of iNeill et al.l (120061 ). 



In the near future, we will use SDSS-II Supernova Survey data to extend this study in 
several directions, including a higher-statistics measurement of the low-redshift rate, mea- 
surement of the SN la rate vs. host galaxy star-formation rate and other host galaxy prop- 
erties, and measurement of the SN la rate to 2; ~ 0.3. The Fall 2006 and Fall 2007 observing 
seasons each yielded ~ 30 spectroscopically confirmed SNe la at redshift z < 0.12, compara- 
ble to the size of the 2005 data set analyzed here. The final SDSS-II Supernova Survey sample 
includes of order 500 spectroscopically confirmed SNe la to 2; < 0.4. 
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